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We update the standard model predictions for the branching ratios of the decays B — » 
(K, K* , X s ) £ + £~, with I = e, fi. Using the measurements of B(B — > X s j) and £?(B — > 
K£ + £~) as well as the upper limits on B(B — * (K* , X s )£ + £~), we work out model indepen- 
dent bounds on the relevant Wilson coefficients. We show the impact of such bounds on the 
parameter space of the MSSM with minimal flavour violation and of supersymmetric models 
with additional sources of flavour changing. 



1 Introduction 

With increased statistical power of experiments at the 5-factories, the exclusive and inclusive 
B — » (K, K* ,X S ) £ + £~ decays will be measured very precisely, in the next several years. On the 
theoretical side, partial results in next-to-next-to-leading logarithmic accuracy are now available 
in the inclusive channels B — > X s £ + £~ 0H. In this talk, we review the status of the standard 
model predictions for the inclusive and exclusive modes and present the model independent 
constraints implied by the new data (see Ref.u for a more detailed presentation of the results). 
The experimental input that we use in our analysis is given below. Except for the inclusive 
branching ratio for B — ► X s ^, which is the average of the results from CLEOQ, ALEPHEl and 
BELLE Q measurements, all other entries are taken from the two BELLE papers listed in Ref.@: 
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B(B -> X s n + H~) < 19.1 x 10" 6 at 90% C.L. , (7) 
B(B -> X s e + e-) < 10.1 x 10~ 6 at 90% C.L. . (8) 

The upper bounds given in Eqs. (p|) - (H) refer to the so-called non-resonant branching ratios 
integrated over the entire dilepton invariant mass spectrum. In the experimental analyses, ju- 
dicious cuts are used to remove the dominant resonant contributions arising from the decays 
B — > (X s , K, K*)( J/ifi, ip' , ...) — > (X s , K, K*)£ + £~. A direct comparison of experiment and 
theory is, of course, very desirable, but we do not have access to this restricted experimental 
information. Instead, we compare the theoretical predictions with data which has been corrected 
for the experimental acceptance using SM-based theoretical distributions. In the present analy- 
sis, we are assuming that the acceptance corrections have been adequately incorporated in the 
experimental analysis in providing the branching ratios and upper limits listed above. We will 
give the theoretical branching ratios integrated over all dilepton invariant masses to compare 
with these numbers. However, for future analyses, we emphasize the dilepton invariant mass 
distribution in the low-s region, s = / rrz^ po ; e < 0.25, where the NNLO calculations for the 

inclusive decays are known, and resonant effects due to J/ip, ip' , etc. are expected to be small. 
The effective Hamiltonian in the SM inducing the b — ► s£ + £~ and b — > S7 transitions is: 

AC 10 

W eff = -T^E C >W°iW ' ( 9 ) 

where Oj(/u) are dimension-six operators at the scale Cj(ju) are the corresponding Wilson 
coefficients, Gf is the Fermi coupling constant, and the CKM dependence has been made explicit. 
The operators can be chosen as in Ref. a 
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where the subscripts L and R refer to left- and right- handed components of the fermion fields. 



2 Expectations in the standard model 

The SM predictions for semileptonic rare B decays are summarized in table |]. 

Note that a genuine NNLO calculation of the inclusive branching ratios only exists for 
values of s below 0.25. For high values of s, an estimate of the NNLO result is obtained 
by an extrapolation procedure. It is possible to show 0, that the full NNLO invariant mass 
distribution is very well approximated, in the entire low-s range, by the partial NNLOj] for the 
choice of the scale /i = 2.5 GeV. This is yet another illustration of the situation often met in 
perturbation theory that a judicious choice of the scale reduces the higher order corrections. It 
seems, therefore, reasonable to use the partial NNLO curve corresponding to /if, = 2.5 GeV as 
an estimate for the central value of the full NNLO for s > 0.25. 

For what concerns the exclusive decays B — * K^£ + £~ , we implement the NNLO corrections 
calculated by Bobeth et al. □ and by Asatrian et al. □ for the short-distance contribution. Then, 



a By partial NNLO we mean that all the terms, for whom the low-s assumption was computationally essential, 
have to be dropped (see Ref. □ for more details) . 



Table 1: SM predictions at NNLO accuracy for the various inclusive and exclusive decays involving the 
quark transition b — > s£ + £~. For the exclusive channels the indicated errors correspond to variations of 
the form factors, /j,b, mt, P oie and m c /mb, respectively. For the inclusive channels the errors correspond, re- 
spectively, to variations of fib, mt.poie and m c /rrib. In the last column we add the errors in quadrature. 



B -► Ki + £~ 


(0.35 ± 0.11 ± 0.04 ± 0.02 ± 0.0005) x 10" 6 


(0.35 ±0.12) x 10" 6 


B -► K*e + e~ 


(l.58 ± 0.47 ± 0.12±g;g| ± 0.04) x 10~ 6 


(1.58 ±0.49) x 10~ 6 


B -► K*fi + fi- 


(l.l9 ± 0.36 ± 0.12t{];{jg ± 0.04) x 10" 6 


(1.19 ±0.39) x 10~ 6 


B -> X sf i+fi- 


(4.15 ± 0.27 ± 0.21 ± 0.62) x 10" 6 


(4.15 ± 0.70) x 10~ 6 


B -» X s e + e~ 


(6.89 ± 0.37 ± 0.25 ± 0.91) x 10" 6 


(6.89 ± 1.01) x 10~ 6 



we use the form factors calculated with the help of the QCD sum rules in Ref. H. In order to 
accommodate present data on the B — > K*j decay, we use the minimum allowed form factors, 
given in Table 5 of Ref. i, as our default set. In our numerical analysis, we add a flat ±15% 
error as residual uncertainty on the form factors. 

Let us stress once more that, from a theoretical point of view, would be much better to 
compare the non-resonant branching ratios integrated in the low-s region with the corresponding 
experimental bounds. Defining such region according to the Belle analysis presented in Ref. 13, 
we choose the integration limits as follows: 



B X,Se- : (°-^Y < i < ( Mj /*-°- 6GeV V . (10) 
\ m b J \ m b J 

B X^V : f < s < ( M ^- a35GeV V . (11) 

For the integrated branching ratios in the SM we find: 

B(B -> X s e + e~) = (2.47 ± 0.40) x 10~ 6 {5B Xa ee = ±16%) , (12) 

low 

B(B -> X^+pT) = (2.75 ± 0.45) x 10~ 6 (5B Xs ^ = ±16%) . (13) 

low 

3 Model independent analysis 

In our analysis, we make the assumption that the dominant new physics effects can be imple- 
mented by using the SM operator basis for the effective hamiltonian. The Wilson coefficients 
that are constrained by the set of measurements given in Eqs. (|l|)-@ are 67,8,9,10- Note that 
b — ► s£ + £~ transitions depend on C7 and Cg only through the effective coefficient A^(Cy, C§)i. 
Our first step consists in writing B(B — » X s j) as a function of Rig^i) = Alfg(fj,)/Ay?$(fj,) 
(following Ref. □) and plotting the regions in the [R7,Rg] plane that are allowed by Eq. ([[]). 
Requiring \Rg(fj,\y)\ — 10 111 order to satisfy the constraints from the decays B — ► X s g and 
B — * Xa (where X^ denotes any hadronic charmless final state) we are, then, able to extract 
the bounds on A-j{^ b ). It was recently pointed out in Ref. Ii9 that the charm mass depen- 
dence of the B — > X s ^f branching ratio was underestimate in all the previous analyses. Indeed, 
the replacement of the pole mass (ni CiPO i e /m by p i e = 0.29 ± 0.02) with the MS running mass 
(m^ 1 (/i b ) /m b}PO i e = 0.22 ± 0.04) increases the branching ratio of about 11%. In order to take 
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Figure 1: 90% C.L. bounds in the [i?7(pi), J?g(/i)] plane following from the world average B —* X a ~/ branching ratio 
for /i = mw (left-hand plot) and \x = 2.5 GeV (right-hand plot). Theoretical uncertainties are taken into account. 
The solid and dashed lines correspond to the m c — m CtPO i e and m c — m* /s (/if,) cases respectively. The scatter 
points correspond to the expectation in MFV models (the ranges of the SUSY parameters are specified in the text). 

into account this additional source of uncertainty, we work out the constraints on the Wilson 
coefficients for both choices of the charm mass; we will then use the loosest bounds in the 
b —* s£ + l~~ analysis. We present the resulting allowed regions in Figs. [l|a and |l|b (note that 
according to the previous discussion we are interested in the bounds at \Xb = 2.5 GeV). The 
regions in Fig. |l|b translate in the following allowed constraints: 

f m c /m b = 0.29 : A 7 ot (2.5 GeV) e [-0.37, -0.18] & [0.21, 0.40] , , , 

{ m c /m b = 0.22 : A 7 ot (2.5 GeV) e [-0.35, -0.17] & [0.25,0.43] . ^ ' 

In the subsequent numerical analysis we impose the union of the above allowed ranges 

- 0.37 < A t 7 ot,<0 (2.5 GeV) < -0.17 & 0.21 < ^ t 7 ° t,>0 (2.5 GeV) < 0.43 (15) 

calling them ^4 7 ot -positive and ^4 7 ot -negative solutions. 

The bounds in the [C^ p {^w)^ Cio P ] plane, implied by the experimental results are 
summarized in Fig. |2| where the two plots correspond respectively to the j4 7 ot -positive and 
A 7 ot -negative solutions just discussed. Note that the overall allowed region is driven by the 
constraints emanating from the decays B — > X s e + e~ (outer contours) and B — ► Kfi + [i~ (inner 
contours). 

4 Expectation in supersymmetry 

In this section we analyze the impact of the b — > 57 and b — ► s£ + l~ experimental constraints 
on several super symmetric models. We will first discuss the more restricted framework of the 
minimal flavour violating MSSM, and then extend the analysis to more general models in which 
new SUSY flavour changing couplings are allowed. 

4-1 Minimal flavour violation 

As already known from the existing literature (see for instance Ref.0), minimal flavour violating 
(MFV) contributions are generally too small to produce sizable effects on the Wilson coefficients 
Cg and Cio- In the MFV scheme all the genuine new sources of flavour changing transitions 
other than the CKM matrix are switched off, and the low energy theory depends only on the 
following parameters: /j,, M2, tan/3, M H ±, and 9^ (see Ref . H for a precise definition of 
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Figure 2: NNLO Case. Superposition of all the constraints. The plots correspond to the A\ ot (2.5 GeV) < and 
A? ot (2.5 GeV) > case, respectively. The points are obtained by means of a scanning over the EMFV parameter- 
space and requiring the experimental bound from B — » Xsj to be satisfied. 



the various quantities). Scanning over this parameter space and taking into account the lower 
bounds on the sparticle masses as well as the b — ► 57 constraint given in Eq. (||), we derive 
the ranges for the new physics contributions to C9 and C\q. In order to produce bounds that 
can be compared with the model independent allowed regions plotted in Fig. ^, we divided the 
surviving SUSY points in two sets, according to the sign of Aj ot : 

A tot < . f qf FV W) G [-0.2,0.4] , 
7 <0 ^lC^G [-1.0,0.7] • (16) 

(Ci IFV (fx w ) G [-0.2,0.3], 
7 I C^ FV E [-0.8,0.5]. (17) 

We stress that the above discussion applies to any supersymmetric model with flavour universal 
soft-breaking terms, such as minimal supergravity MSSM and gauge-mediated supersymmetry 
breaking models. Beyond-the-SM flavour violations in such models are induced only via renor- 
malization group running, and are tiny. Hence, they can be described by MFV models discussed 
in this paper. 

Before finishing this subsection and starting our discussion on models with new flavour 
changing interactions, let us show in more detail the impact of b — > sj on MFV models. The 
scatter plot presented in Fig. |l| is obtained varying the various MFV SUSY parameters and 
shows the strong correlation between the values of the Wilson coefficients C7 and Cs- In fact, 
the SUSY contributions to the magnetic and chromo-magnetic coefficients differ only because 
of colour factors and loop-functions. 



4-2 Chargino contributions: Extended-MFV models 

A basically different scenario arises if chargino-mediated penguin and box diagrams are con- 
sidered in connection with non-vanishing mass insertions (See Ref. for a definition of the 
so-called mass insertion approximation (MIA)). As can be inferred by Table 4 in Ref. 1111, the 
presence of a light ti generally gives rise to large contributions to C9 and especially to C\q. In 
the following, we will concentrate on the so-called Extended MFV (EMFV) models introduced 
in Ref. Ej. In these models we can fully exploit the impact of chargino penguins with a light t 
still working with a limited number of free parameters. EMFV models are based on the heavy 
squarks and gluino assumption. In this framework, the charged Higgs and the lightest chargino 
and stop masses are allowed to be light while the rest of the SUSY spectrum is assumed to be 



degenerate and heavier than 1 TeV. The assumption of a heavy gluino suppresses any possible 
gluino-mediated SUSY contribution to low energy observables. Note that even in the presence of 
a light gluino these penguin diagrams remain suppressed due to the heavy down squarks present 
in the loop. In the MIA approach, a diagram can contribute sizably only if the inserted mass 
insertions involve the light stop. This leaves us with only two unsuppressed flavour changing 
sources other than the CKM matrix, namely the mixings ul — ti (denoted by 5u L t 2 ) a f4 ~ ^2 
(denoted by &c L t 2 )- ^he phenomenological impact of ^i^ L ^ as been studied in Ref. EJ and its 
impact on the b — > 57 and b — > s£ + £~ transitions is indeed negligible. Therefore, we are left 
with the MIA parameter <5f 2 g, only. Scanning over the SUSY parameter space (/U, M2, tan/3, 
Mj 2 , sin#£, Mg±, My and <5t 2 c L ) an d imposing the constraints from the sparticle masses lower 
bounds and B — > X s ^, we obtain the points plotted in Fig. Note that these SUSY models 
can account only for a small part of the region allowed by the model independent analysis of 
current data. 



5 Conclusions 



We have presented theoretical branching ratios, model independent analyses and SUSY pre- 
dictions for the rare B decays B — > X s £ + £~ and B — > (K, K*)£ + £~ , incorporating the NNLO 
improvements. The dilepton invariant mass spectrum is calculated in the NNLO precision in the 
low dilepton invariant mass region, s < 0.25 and an extrapolation is used for s > 0.25. Current 
B factory experiments will soon probe these decays at the level of the SM sensitivity and we 
stress the need to measure the inclusive decays B — > X s £ + £~ in the low dilepton mass range for 
a proper comparison with the SM expectations. 
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